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We discuss the exact solution of Einstein's equation corresponding to a static and plane 
symmetric distribution of matter with constant positive density located below z = 
matched to vacuum solutions. The internal solution depends essentially on two constants: 
the density p and a parameter re. We show that these space-times finish down below at 
an inner singularity at finite depth d < \p^ > ■ We show that for re > 0.3513 . . . , the 
dominant energy condition is satisfied all over the space-time. 

We match these singular solutions to the vacuum one and compute the external 
gravitational field in terms of slab's parameters. Depending on the value of re, these 
slabs are either attractive, repulsive or neutral. The external solution turns out to be a 
Rindler's space-time. Repulsive slabs explicitly show how negative, but finite pressure 
can dominate the attraction of the matter. In this case, the presence of horizons in the 
vacuum shows that there are null geodesies which never reach the surface of the slab. 

We also consider a static and plane symmetric non-singular distribution of matter 
with constant positive density p and thickness d (0 < d < surrounded by two 

external vacuums. We explicitly write down the pressure and the external gravitational 
fields in terms of p and d. The solution turns out to be attractive and remarkably asym- 
metric: the "upper" solution is Rindler's vacuum, whereas the "lower" one is the singular 
part of Taub's plane symmetric solution. Inside the slab, the pressure is positive and 
bounded, presenting a maximum at an asymmetrical position between the boundaries. 
We show that if < ^/Wirpd < 1.527 . . . , the dominant energy condition is satisfied all 
over the space-time. We also show how the mirror symmetry is restored at the Newtonian 
limit. 

We also find thinner repulsive slabs by matching a singular slice of the inner solution 
to the vacuum. 

We also discuss solutions in which an attractive slab and a repulsive one, and two 
neutral ones are joined. We also discuss how to assemble a "gravitational capacitor" by 
inserting a slice of vacuum between two such slabs. 



1. Introduction 

Due to the complexity of Einstein's field equations, one cannot find exact solutions 
except in spaces of rather high symmetry, but very often with no direct physical 
application. Nevertheless, exact solutions can give an idea of the qualitative features 
that could arise in General Relativity, and so, of possible properties of realistic 
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2 Ricardo E. Gamboa Saravi 

solutions of the field equations. 

We have recently discussed exact solutions of E inste in's equation presenting an 
empty (free of matter) singular repelling boundary \U21. These singularities are not 
the sources of the fields, but they arise owing to the attraction of distant matter. 

In this paper, we want to illustrate this and other curious features of relativistic 
gravitation by means of a simple exact solution: the gravitational field of a static 
plane symmetric relativistic perfect incompressible fluid with positive density lo- 
cated below z — matched to vacuum solutions. In reference^, we analyze in detail 
the p roperties of this internal solution, originally found by A. H. Taub SI (see also 
| 5 | 6 | 7 |^ and we find that it finishes up down below at an inner singularity at finite 




depth d, where < d < y 2ip" Depending on the value of a parameter k, it turns out 
to be gravitational attractive (k < K cr it), neutral (k = Kent) or repulsive (re > K C rit), 
where K cr it = 1.2143 .... We also show that for re > 0.3513 . . . , the dominant energy 
condition is satisfied all over the space-time. 

In this paper, we make a detailed analysis of the matching of these exact solutions 
to vacuum ones. Here, we impose the continuity of the metric components and of 
their first derivatives at the matching surfaces, in contrast to reference^, where not 
all these derivatives are continuous at the boundary. 

In the first place, we consider the matching of the whole singular slabs to the 
vacuum, and explicitly compute the external gravitational fields in terms of the 
slab parameters. Repulsive slabs explicitly show how negative but finite pressure 
can dominate the attraction of the matter. In this case, they have the maximum 
depth, i.e., d = y 1 '^Ip' anc ^ ex t er ior solution presents horizons showing that 
there are vertical photons that cannot reach the slab surface. 

Secondly, we consider a non-singular slice of these slabs with thickness d (0 < 
d < . Tjf - ) surrounded by two external vacuum. Some of the properties of this 



solution have already been discussed in reference^. Here, we explicitly write down 
the pressure and the external gravitational fields in terms of p and d. The solution 
turns out to be attractive, and remarkably asymmetric: the "upper" solution is 
Rindler's vacuum, whereas the "lower" one is the singular part of Taub's plane 
symmetric solution. Inside the slab, the pressure is positive and bounded, presenting 
a maximum at an asymmetrical position between the boundaries. We show that if 
< ^/6npd < 1.527..., the dominant energy condition is satisfied all over the 
space-time. This solution finishes up down below at an empty repelling boundary 
where space-time curvature diverges. This exact solution clearly shows how the 
attraction of distant matter can shrink the space-time in such a way that it finishes 
at a free of matter singular boundary, as pointed out in M \y e also show how the 
mirror symmetry is restored at the Newtonian limit. 

We also construct thinner repulsive slabs by matching a singular slice of the 
inner solution to vacuum. These slabs turn out to be less repulsive than the ones 
discussed above, since all incoming vertical null geodesies reach the slab surface in 
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The Gravitational Field of a Plane Slab 3 

this case. 

For the sake of completeness, in section[2] we include some results from reference 
^ which are necessary for the computations of the following sections. In section[3j we 
show under which conditions the dominant energy condition is satisfied. In section 
|4l we discuss how solutions can be matched. In section [5l we study the matching of 
the whole singular interior solution to vacuum. In section [6l we match two interior 
solutions facing each other. In section [71 we discuss the matching of a non-singular 
slice of the interior solution with two different vacuums. In section [51 we show how 
the mirror symmetry of this solution is restored at the Newtonian limit. In section 
[9] we construct thinner repulsive slabs [d < y^^) by matching a singular slice of 

the inner solution to vacuum. 

Throughout this paper, we adopt the convention in which the space-time metric 
has signature (— + H — h), the system of units in which the speed of light c = 1, 
Newton's gravitational constant G = 1 and g denotes gravitational field and not 
the determinant of the metric. 



2. The interior solution 

In this section, we consider the solution of Einstein's equation corresponding to a 
static and plane symmetric distribution of matter with constant positive density 
and plane symmetry. That is, it must be invariant under translations in the plane 
and under rotations around its normal. The matter we shall consider is a perfect 
fluid of uniform density p. The stress-energy tensor is 

Tab = {p + p) u a u b + p g ab , (1) 

where u a is the velocity of fluid elements. 

Due to the plane symmetry and staticity, following ® we can find coordinates 
(t, x, y, z) such that 

ds 2 = -Q{zf dt 2 + e 2V ^ (dx 2 + dy 2 ) + dz 2 , (2) 

that is, the more general metric admitting the Killing vectors d x , d y , xd y — yd x and 
d t . 

The non identically vanishing components of the Einstein tensor are 

G tt = -Q 2 {2V" + 3V' 2 ) (3) 
G xx = G yy = e 2V (g"/Q + Q'/Q V + V" + V 12 ) , (4) 

G zz = V (2 Q'/Q + V) , (5) 

where a prime (') denotes differentiation with respect to z. 

On the other hand, due to the assumed symmetries and to the fact that the 
material content is a perfect fluid, u a = {—Q, 0, 0, 0), so 



T ab = dia,g(pQ 2 ,pe 2V ,pe 2V ,p) , 



(0) 
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4 Ricardo E. Gamboa Saravi 

where p depends only on the z-coordinante. Thus, Einstein's equations, i.e., G ab = 
8iTT ab , are 

2V" + 3V' 2 = -87rp, (7) 
Q"/Q + G'/G V + V" + V' 2 = 8tt P , (8) 
V'(2g'/g + V')=8Trp. (9) 

Moreover, V ' a T ab = yields 

p' = -{ P + p)G'/G. (10) 

Of course, due to Bianchi's identities equations, (0), j8]), j9|) and (fTUf are not inde- 
pendent, so we shall here use only (jTj) , (J9J) , and (fT0|) . 
Since p is constant, from (|10|) we readily find 

P = c p /g(z)-p, (ii) 

where C p is an arbitrary constant. 

By setting W(z) = e^ v< ^ z \ we can write ([7]) as W" = —6npW, and its general 
solution can be written as 

W(z) = Cx sin (yfthrp z + C 2 ), (12) 

where C\ and C 2 are arbitrary constants. Therefore, we have 

V{z) = ^ In (Ci sin (y/&irp z + C 2 )) . (13) 

Now, by replacing (fTT) into ((9]), we get the first order linear differential equation 
which G{z) obeys 

<r~(£ + £)« + ^ 

= — \J 6np [tan u -\ — cot u ) G + \J 67rp — tan it, (15) 
V 3 / P 

where u — y/6irp z + C 2 - And in the last step, we have made use of (fT3|) . The general 
solution of ([14]) can be written as @ 

(smu) i / J \ P Jo (cosu') 2 / 

_ cosu 3C„ . 2 _/ 2 13 . o \ 

= C 3 - + — - sm 2 u 2 F 1 1,-;— ;sm 2 w , 16 

(sinw)^ 7 P V 3 6 / 

where C3 is another arbitrary constant, and 2-P1 (a, b; c; 0) is the Gauss hypergeo- 
metric function (see the appendix at the end of the paper). 
Therefore, the line element <j2j) becomes 

ds 2 = -G{z) 2 dt 2 + (d smu)^ (dx 2 + dy 2 ) + dz 2 , (17) 



a In the appendix we show how the integral appearing in the first line of l|16| l is performed. 
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The Gravitational Field of a Plane Slab 5 

where Q(z) is given in (|16p and u = \/Qnp z + C%. Thus, the solution contains five 
arbitrary constants: p, C p , C\, C2, and C3. The range of the coordinate z depends 
on the value of these constants. 

Notice that the metric (|17[) has a space-time curvature singularity where sin u = 
0, since straightforward computation of the scalar quadratic in the Riemann tensor 
yields 

R abcd R abcd = 4 (g" 2 + 2 g' 2 v' 2 ) /g 2 + a (2 v" 2 + 4 v"v' 2 + 3 v' A ) 

256 2 2 /o , „• -4.. , 3 /p , ,\/ 3 ? 



3 *n 2+8i - "Mrvlf-vJ' <18) 

so R a bcdR ahcd — * 00 when sinw — > 0. 

On the other hand, by contracting Einstein's equation, we get 

R{z) = 8tt(p - 3p(z)) = 87r(4p - 3C P /G{z)) . (19) 

For p > 0, C p > and C 3 > 0, the solution ((TTJ) was found by Taub EE. 
Nevertheless, this solution has a wider range of validity. 

Of course, from this solution we can obtain vacuum ones as a limit. In fact, when 
C p = 0, it is clear from (fTTj) that p(z) = —p, and the sol ution fj 1 7|) turns out to be 
a vacuum solution with a cosmological constant A = 87rp 

ds 2 — —cos 2 u sin~^ udt 2 + sin 5 u (dx 2 + dy 2 ) + dz 2 , 
— oo < t < oo, — oo < x < oo, — oo < y < oo, < u < 7r, (20) 

where u = A/3A/2 z + Ci. We get from p9|) that it is a space-time with constant 
scalar curvature 4A, and from (| 18[) we get that 



ibcd ^ 2 ' ~ ' 



RabcdR = - AM 2 + — j- . (21) 

o \ sin 11/ 

Now, we take the limit A — > (p — > 0). By setting C2 = 7r — -^^p and an 
appropriate rescaling of the coordinates {£, x, y}, we can readily see that, when 
A->0. ([20)1 becomes 

ds 2 = -(1 - 3gz)~i dt 2 + (1 - 3.gz)s (dx 2 + cfy 2 ) + dz 2 , 
—00 < t < 00, — 00 < a; < 00, — 00 < y < 00, < 1 — 3gz < 00 , (22) 

where g is an arbitrary constant. In (|22p . the coordinates have been chosen in such 
a way that it describes a homogeneous gravitational field g pointing in the negative 
z-direction in a neighborhood of z = 0. The metric ([22)) is Taubs's ^ vacuum plane 
solution expressed in the coordinates used in Ref. where a detailed study of it 
can be found. 

On the other hand, by setting C2 = § + and an appropriate rescaling of 
the coordinate t, we can readily see that, when A — > 0, (|20"|) becomes 



ds 2 = -(1 + gz) 2 dt 2 + dx 2 + dy 2 + dz 2 , 

1 , . 

-00 < t < 00, —00 < x < 00, —00 < y < 00, < z < 00 , (23) 

9 
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where g is an arbitrary constant, and the coordinates have been chosen in such a way 
that it also describes a homogeneous gravitational field g pointing in the negative 
z-direction in a neighborhood of z = 0. The metric (|23l) is, of course, Rindler's fiat 
space-time. 

For exotic matter, some interesting solutions also arise, but the complete anal- 
ysis turns out to be somehow involved. So, for the sake of clarity, we shall confine 
our attention to positiv e va lues of p and C p ^ 0, leaving the complete study to a 
forthcoming publication HI. 

Now, it is clear from ([7]), (jHJ), and (flT)]) that field equations are invariant under 
the transformation z — > ±z + zq, i.e., z-translations and mirror reflections across 
any plane z =const. Thus, if {Q[z), V(z),p(z)} is a solution {G(±z + z ), V(±z + 
zq),p(±z + zq)} is another one, where zq is an arbitrary constant. Therefore, taking 
into account that u = \/6np z + C2 , without loss of generality, the consideration of 
the case < u < tt/2 shall suffice. 

By an appropriate rescaling of the coordinates {x, y}, without loss of generality, 
we can write the metric (fTTj) as 

ds 2 = -Q{zf dt 2 + sin^ u (dx 2 + dy 2 ) + dz 2 , 
—00 < t < 00, —00 < x < 00, —00 < y < 00, < u = z + C2 < 7r/2, 

(24) 



and (fT6j) as 

nC p cosu 3C„ . 2 / 2 13 . 2 \ 

G(z) = — p - — + —P. S m 2 u 2F1 1,-;— ;sm 2 w , (25) 

P sins u "P V 3 6 I 

where k is an arbitrary constant. 

By replacing (|2"5j) into (jTTJ) , we see that the pressure is independent of C p . On 
the other hand, since Q(z) appears squared in (f24|) . it suffices to consider C p > 0. 
Therefore, rescaling the coordinate t, we may set C p = p. Thus, (f2"5|) becomes 

^/ s ^ , •, cosu 3., „ /„ 2 13 . 9 \ . 

g z = G K (u) =k + - sm 2 M 2 Fi 1,-;— ;sm 2 W , 26 

sin 1/J u 7 V 3 6 / 

where G K (u) is defined for future use, and we recall that u = yf&tzp z + C2. Further- 
more, (jllj) becomes 

p(z) = />(l/a(z)-l). (27) 

Therefore, the solution depends on two essential parameters, p and K. We shall 
discuss in detail the properties of the functions Q(z) and p(z) depending on the 
value of the constant k. 

By using the transformation (|93j) . we can write G(z) as 

(k - K cx it) cosu + 2 F 1 (-^, -i; i;cos 2 u ) 
5(z) = rrji L , (28) 
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where 

\/7rr(7/6) , , 

^ = ^^ = 1.2143..., (29) 

which is the form used in references and which is more suitable to analyze its 
properties near u = tt/2. 

Now, the hypergeometric function in (|26j) is a monotonically increasing contin- 
uous positive function of u for < u < 7r/2, since c — a — 6 = 1/2 > 0. Furthermore, 
taking into account that 2-Fi(a, b; c; 0) = 1 and (|93|) . we have 

Therefore, we readily see from (|26[) that, no matter what the value of k is, 
G{ z )\u=ir/2 — 1) an d we get then from (|27|) that p(z) vanishes at u = ir/2. On 
the other hand, since 

5(z) = K!i^ + 0(J) as u->0, (31) 

£(z)U=o = if k = 0, whereas it diverges if ft ^ 0. 

For the sake of clarity, we shall analyze separately the cases k > 0, k — 0, and 

K < 0. 



2.1. /« > 

In this case, it is clear from (|26|) that C/(z) is positive definite when < u < ir/2. 
On the other hand, from |(5J) and ([9]), we get 

S" = g'v'~ gv" = - (v" + ^- + 4^) g + AttC p = v' 2 g + 4tt P , (32) 

where we have made use of (IT41 . ([7]) and C p = p. Then, also g" is positive definite 
in < u < 7r/2, and so g' is a monotonically increasing continuous function of u in 
this interval. 

Now, taking into account that g' — d z g — y/WrTp d u g, we get from ([26| that 

g'( z ) = _^v|^ u -| + 0( u |) asu ^o, (33) 

and from (gSJ) that 

0'(2!)U=7r/2 = V 6?r P ( K crit - «) • (34) 

If k > K cr it, 5' is negative for small enough values of u and non-positive at 
u = n/2. Hence g 1 is negative in < u < tt/2, so g(z) is decreasing, and then 
g(z) > g(z)\ u=7T /2 = 1 in this interval (see FigfjTa) an d FigQJb)). 

For K cr it > k > 0, g 1 is negative for sufficiently small values of u and positive at 
tt/2. So, there is one (and only one) value u rn where it vanishes. Clearly g(z) attains 
a local minimum there. Hence, there is one (and only one) value uo (0 < uq < tt/2) 
such that g{z)\ u=Uo — g(z)\ u=7T /2 = 1? an d then g(z) < 1 when uo < u < tt/2 (see 
FigHlJc) and FiglUd)). 
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8 Ricardo E. Gamboa Saravi 





Since Q(z) > 0, it is clear from (|2"T|) that p(z) > if Q{z) < 1, and p(z) reaches 
a maximum when Q{z) attains a minimum. 

Therefore, for k > K CI it, p(z) is negative when < u < n/2 and it increases 
monotonically from — p to and it satisfies \p\ < p all over the space-time (see 
FigOJa) and Figrjjb)). 

On the other hand, for K cr it > n > 0, p(z) grows from — p to a maximum positive 
value when u = u m where it starts to decrease and vanishes at u = Tt/2. Thus, p{z) 
is negative when < u < uq and positive when uq < u < tt/2 (see Figfljc) and 
FigfJJd)). It can be readily seen from PS)) and (|27[) that, as k decreases from K cr i t to 
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1 




tt/4 tt/2 tt/4 tt/2 

(a) re = (b) re < 



Fig. 2. §(2) and p(z) as functions of u for k < 0. 

0, u m moves to the left and the maximum value oip(z)/p monotonically increases 
from to 00. In section [3J we shall show that for re = Kd ec — 0.3513 ... it gets 1, 
and then for < re < Kdec, there is a region of space-time where p > p and where 
the dominant energy condition is thus violated. 

2.2. k — 

In this case, it is clear from (1261) that Q monotonically increases with u from 
to G{z)\ u=7! /2 = 1- Therefore, p is a monotonically decreasing positive continuous 
function of u in < u < tt/2 (see FiglJJa)). Furthermore, at u = it diverges, since 

p{z) ~ -7r u ~ 2 ~ * +°° as u ~~ * 0. (35) 

2.3. k < 

In this case, we see from (f3"3"l) that 5' is positive when u takes small enough values, 
and from (|34|) we see that it is also positive when u is near to tt/2. 

Now, suppose that Q' (z) attains a local minimum when ti = U\ (0 < Ui < tt/2), 
then 5"(z)| u=Ul = 0. Hence, we get from (|32|) that G(z)\ u=Ul < 0. And taking into 
account that V'(z)\ u=Ul = 2y/6irp/3 cotui > 0, we see from (TTJ} that 5'(z)| M=Ml > 
0. Thus, we have shown that G'(z) is a continuous positive definite function when 
< u < tt/2 if k < 0. 

Therefore, in this case, £7(z) is a continuous function monotonically increasing 
with u when < u < tt/2. Since it is negative for sufficiently small values of u and 1 
when u = tt/2 it must vanish at a unique value of z when u — u K (say). Furthermore 
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G(z) < 1 when < u < ir/2. Clearly, we get from (|2"8|) that u K is given implicitly in 
terms of k through 

K = Kcrit 77^ • (36) 

cosm k sin ' u K 

We can readily see from (|36|) that u K is a monotonically decreasing function of k in 
— oo < k < 0, and it tends to ir/2 when k — > — oo and to when k — ► _ . 

From (|27|) . it is clear that p(z) diverges when u = u K . Furthermore, ([27]) also 
shows that p(z) < when G(z) < 0. And taking into account that G(z) < 1, we see 
that p{z) > when G{z) > 0. Therefore, p{z) is negative when < u < u K whereas 
it is positive when u K < u < ir/2 (see Fig(2][b)). 

On the other hand, we see from (fT8|) that, when k is negative, another space- 
time curvature singularity arises at u K (besides the one at u = 0) since p diverges 
there. 

Therefore, if k is negative, the metric (|24|) describes two very different space- 
times: 

(a) For < u < u K , the whole space-time is trapped between two singularities 
separated by a finite distance ^/Qi:pu K . This is a space-time full of a fluid with 
constant positive density p and negative pressure p monotonically decreasing with 
u, and p(z)\ u= o = —p and p(z) — > — oo as u — * u K . 

(b) For u K < u < 7r/2, the pressure is positive and monotonically decreasing 
with u, p(z) — > oo as u — * m k and p(2:)| u=7 r/2 = 0. 

3. The maximum of the pressure and the dominant energy 
condition 

We have seen in the preceding section that for k > K CI it, p(z) is negative, it increases 
monotonically from — p to and it satisfies \p\ < p all over the space-time (see 
FigHIa) and Fig[T](b)). On the other hand, for k < 0, p(z) is unbounded at an inner 
singularity and thus the dominant energy condition is not satisfied in this case. 

For K clit > k > 0, since Q(z) > 0, it is clear from (fTTj) that p{z) > if G{z) < 1, 
and that p(z) reaches a maximum when G{z) attains a minimum. Then, p(z) grows 
from — p to a maximum positive value p m when u = u m , where it starts to decrease 
and vanishes at u — n/2. Thus, —p < p(z) < for < u < uq and < p(z) < p m 
when u n < u < ir/2 (see FigQJc)). 

We readily see from ((9]), since G'(z)\ u=Um vanishes, that 

Pm = p(z)\ u =u m = — {V'{z)f \ u =u m = ~ cot 2 u m , (37) 

87T 3 

where we have made use of (|13[) , and so the maximum value of p(z) monotonically 
decreases from oo to in < u rn < tt/2. 

Now, by replacing (|37| into (]28[) and taking into account (|27ll . we can write 
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12 3 

Fie. 3. K ELS cL function of p m . 



down k in terms of p m 

(3p m )^ V(p + Pm)(p + 3p ro )B 2 2 ' 6 1 2 1 p + 3 Pm )j 

_ (38) 

o fP™ ( \ 1 Pm . 8 p m Z \ 
= ^-2^—[l--— + --^ + ...j iorp m <p, 

(39) 

which clearly shows that n — > k ct h as p m — > 0. On the other hand, by using (|96[) . 
we can write 

which clearly shows that k — > as p m — > oo. Thus, as k increases from to K C rit, 
p m monotonically decreases from oo to (see Figj3|). 

Hence, there is a value Kdoc of k for which p m = p, and from (|38| we see that it 
is given by 

« dec = « cri t + ^ - »*(--, - J; 5; I)) = 0-351307 .... (41) 

Also note that, in this case, we get from ()37|1 that the maximum of the pressure 
occurs at u m = it/6. 

Thus, for < k < ftdco there is a region of space-time where p > p and 
where the dominant energy condition is thus violated. However, we see that for 
Nice < « < Kcrit, the condition \p\ < p is everywhere satisfied. 

Therefore, the dominant energy condition is satisfied all over the space-time if 

K 2^ ^dec- 

Notice that, for K cr it > k > 0, by eliminating k by means of (|38p . the solution 
can be parameterized in terms of p m and p. 
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4. The matching of solutions and the external gravitational fields 

We shall discuss matching the interior solution to a vacuum one, as well as joining 
two interior solutions facing each other at the surfaces where the pressure vanishes. 
For any value of K, p(z) — at u = 7r/2, while for K cr it > k > it also vanishes at 
u = uq. Therefore, the matching at u — 7r/2 is always possible, while the matching 
at u = uq is also possible in the latter case. 

We shall impose the continuity of the metric components and of their first deriva- 
tives at the matching surfaces. 

Notice that, due to the symmetry required, vacuum solutions satisfy the field 
equations (O, ((8]) and ((9]), with p = p = 0. In this case, we immediately get from 
® that or V = or 2 Q'jQ + V' = 0. 

In the former case, we get from ([8|) that Q" = 0, and the solution is 

ds 2 = -(A + Bzf dt 2 + C(dx 2 + dy 2 ) + dz 2 , (42) 

which is the Rindler space-time. 

In the latter one, it can be written as 

ds 2 = -(A + Bz)~i dt 2 + C(A + Bz)% (dx 2 + dy 2 ) + dz 2 , (43) 

which is the Taub's vacuum plane solution^. 

Therefore, as pointed out by the authors of reference ®, if at the matching 
"plane" the interior V' vanishes, we can only match it to Rindler's space-time, 
since for the Taub's one V does not vanish at any finite point. Whereas, if on the 
contrary, V' does not vanish at the matching "plane" , we can only match the inner 
solution with Taub's one, since for Rindler's one, V vanishes anywhere. 

Now, we see from (|T3j) that V' vanishes at u = tt/2 and it is non zero at 
u = uq n/2. Therefore the solution can be matched to Ridler's space-time at 
u = tt/2 and to Taub's vacuum plane solution at u = uq. 

Notice that in reference ^ we did not demand the continuity of V'(z) at the 
matching surface and we analyzed there the matching of the solution to Taub's 
vacuum plane solution at u = tt/2. 

In the next section, we discuss the matching of the whole interior solution to 
Rindler vacuum, for any value of k at u = tt/2, while we match two interior solutions 
facing each other at u = tt/2 in section [6j 

In section[7l for K cr i t > n > 0, we discuss the matching of the slice of the interior 
solution uq < u < tt/2 with both vacua, while, in section [9] we match the remaining 
piece (i.e. < U < uq) to a Taub's vacuum. 

5. Matching the whole slab to a Rindler space-time 

In this section, we discuss matching the whole interior solution to a vacuum one at 
u = tt/2. 

Since the field equations are invariant under z-translation, we can choose to 
match the solutions at z = without losing generality. So we select Ci — tt/2, and 
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then 



becomes 



Q(z) = G K ( x /(^p~z + ir/2) 
- (k - K crit ) sm(y/6irpz) + 2-F1/-5, -f ; |; sin 2 ( % /67rpz)^ 
cos 1 / 3 (v / 67rpz) 

Therefore, the metric (|24)) reads 

ds 2 = -Q(z) 2 dt 2 + cos^ (V67rp z) (dx 2 + dy 2 ) + dz 2 



(44) 



-oo < t < oo, — oo < a; < oo, 



-co < y < oo, 




< z < 0. 



(45) 



We must impose the continuity of the components of the metric at the matching 
boundary. Notice that g tt (0) = -£(0) 2 = -1, g xx {0) = g yy (0) = 1, and p(0) = 0. 

Furthermore, we also impose the continuity of the derivatives of the metric 
components at the boundary. From (|34p . we have 



intcri 

and, from ([45)) we get 



-2 G(0) G'(0) = -2 (Kelt - k) , 



(46) 



9 z g xx (0)| interior = <9zffyy(0)| interior = 5 COS 3 ( y/6np z) SUv(\/ 6n p z) 



The exterior solution, i.e. for z > 0, is the Rindler space-time 

ds 2 = -(1 + gz) 2 dt 2 + dx 2 + dy 2 + dz 2 , 



= 0. 



(47) 



— 00 < t < 00, —00 < x < 00, —00 < y < 00, < z < 00 , (48) 

which describes a homogeneous gravitational field —g in the vertical (i.e., z) direc- 
tion. 

Since g t t(0)\ exterior = -1 and s- xx (0)| exterior = SW (°) I exterior = 1, the continuity 
of the metric components is assured. And, concerning the derivatives, we have 

&z gxx (,z) I exterior — &zgxx (-^) | exterior — 0i (^9) 

which identically matches to 1(17)). 
Moreover, we readily get 

9 z gtt (z)|extcrior = (1 + 2ffz) . (50) 

Then, by comparing it with (133j) , we see that the continuity of d z gu at the boundary 
yields 

g = y/ 6irp (Kcrit - k) , (51) 
which relates the external gravitational field g with matter density p and k. 
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Case 


K 


9 


p(z) 


\P\ <P 


Depth 


Fig. 


I 


K > ftcrit 


< 


-P < p(z) < 


yes 


/ 7T 

V 24p 


Ola) 


II 


K = K C rit 


= 


-p < p{z) < 


yes 


/ TT 


ntb) 


III 


^crit > ft 2^ ^dcc 


> 


-P < P(z) < Pm(«) < P 


yes 


/ 7T 

Y 24p 




IV 


Kdoc > « > 


> 


~P < p(z) < Pm(«) 


no 


/ 7T 

V 24 P 


ntd) 


V 


> K 


> 


unbounded 


no 


(tt/2— ««) 
\/67rp 


El 



Now, by replacing k from (j5Tj) into l]44p. we get 

5 sin(Y / 67rpz) + \f&Tip 2-Fif - 31 — |j sin 2 (-y/ 677/3 2;) 



/67rp cos 1 / 3 (v / 67T/9z) 



(52) 



and the solution is parameterized in terms of the external gravitational field g and 
the density of the matter p. 

It can readily be seen from (|51|) that, if K > K C rit, g is negative and the slab 
turns out to be repulsive. If K — K C rit it is gravitationally neutral, and the exterior 
is one half of Minkowski's space-time. If n < K cr ; t , it is attractive. 

If k > 0, the depth of the slab is ^J^; independently of the value of k. In this 
case, the pressure is finite anywhere, but it is negative deep below and p = —p 
at the inner singularity (see FigQJa), FigQJb) and FigQTc)). But, as discussed in 
section [3J only when k > Kdec is the condition \p\ < p everywhere satisfied. 

If k < 0, the pressure inside the slab is always positive, and it diverges deep 
below at the inner singularity (see FigE]). Its depth is 



d = {tt/2 - u K )/y/6irp, 



(53) 



where u K (0 < u K < tt/2) is given implicitly in terms of k through (|36p . By using 
(|36p . we can write k in terms of d 



' "V 2 ' 6 ' 2 ' * ' •) lr.A\ 

K = Kcrit ~ 1, — ■ (54) 



; sin 2 (^/67rp d) 



sin( v / 67rp d) cos 3 (^JQirp d) 

Now, in this case, by using (|5ip we can write the external gravitational field g in 
terms of the matter density p and the depth of the slab d 

VGtTp _/ 1 11 



sin(v / 67r/5 d) cos 3 (y/6np d) 



iF^--,--; sin 2 (y^pd) 



(55) 



For the sake of clearness, we summarize the properties of the solutions discussed 
above in Table ??. 
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Case 


Quadrant 


T 


Z 


-dT 2 + dZ 2 


Attractive 


I 


(z + 1/g) sinh gt 


(z + 1/g) coshgi 


-{z + l/gfdt 2 + dz 2 


Repulsive 


I and II 


(z-l/g) sinh gt 


(z — 1/g) coshgt 


-(z- l/g) 2 dt 2 +dz 2 


III 


(z — 1/g) coshgt 


(z — 1/g) sinhc/i 


-dz 2 + (z-l/g) 2 dt 2 


IV 


(1/g — z) coshgt 


(z — 1/g) sinhc/i 


-dz 2 + (z-l/g) 2 dt 2 



Some remarks are in order. First, notice that the maximum depth that a slab 
with constant density p can reach is v/^^j being the counterpart of the well-known 

bound M < 4i?/9 (R < ) , which holds for spherical symmetry. 

If we restrict ourselves to non "exotic" matter, the dominant energy condition 
will put aside cases IV and V, as shown in section[3J However, as already mentioned, 
it is satisfied for the cases I, II and III (see Figflja), FigJTJb) and Fig[ljc)). Thus, 
there are still attractive, neutral and repulsive solutions satisfying this condition. 

In this case, we readily get from (|5Tj) the bound 

9 < \J^P (Kcrit - K d ec) » 3.75 y/p . (56) 

In order to analyze the geodesies in the vacuum, it is convenient to consider 
the transformation from Rindler's coordinates t and z to Minkowski's ones T and 
Z shown in Table ??. Notice that, for the repulsive case, four Rindler's patches 
are necessary to cover the whole exterior of the slab. Also note that, in this case, 
z becomes the temporal coordinate in quadrants III and IV, see Fig. 21 In this 
coordinates, of course, the vacuum metric becomes 

ds 2 = -dT 2 + dx 2 + dy 2 + dZ 2 . (57) 

Notice that, the "planes" z = constant correspond to the hyperbolae Z 2 — T 2 — 
constant, and t — constant. On the other hand, incoming vertical null geodesies are 
Z + T = constant, and outgoing ones are given by Z — T = constant. 

For attractive slabs, we readily see from Fig. [Ha) that all incoming vertical 
photons finish at the surface of the slab, while all outgoing ones escape to infinite. 
Vertical time-like geodesies start at the surface of the slab, reach a turning point 
and fall down to the slab in a finite amount of coordinate time t. Notice that a 
particle world-line is tangent to only one hyperbola Z 2 — T 2 = C, with C > 1/g, 
and that the maximum value of z that it reaches is C — 1/g. 

For repulsive slabs, Fig. HJb), two horizons appear in the vacuum: the lines 
T = ±Z, showing that not all the vertical null geodesies reach the surface of the 
slab. In fact, only vertical incoming photons coming from region IV end at the slab 
surface, and only the outgoing ones finishing in region III start at the slab surface. 
Incoming particles can reach the surface or bounce at a turning point before getting 
it. 
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(a) Attractive slab (b) Repulsive slab 



Fig. 4. Vertical time-like and null geodesies in the vacuum 

6. Matching two slabs 

Now we consider two incompressible fluids joined at z = where the pressure 
vanishes, the lower one having density p and the upper having density p' . Thus, the 
lower solution is given by (|45[) . By means of the transformation z — ► —z, p — * p' 
and k — > k' we get the upper one 

ds 2 = -G«/(tt/2- ^Kp'zfdt 2 + cos^(v/f3Vz) (dx 2 + dy 2 ) +dz 2 , 
— oo <t< co, — oo < x < co, — oo < y < co, < z < \j 2^ i ' 

From (I45| . (T47|) and (|58[). we can readily see that gtt{z), g X x{z) and d z g xx (z) are 
continuous at z = 0. Furthermore, from ()46|) we see that the continuity of d z gu 
requires 

y/p (K cr it - K) = -\fp' (K C rit ~ «') ■ (59) 

Thus, if one solution has a k greater than K cr ;t, the other one must have it smaller 
than K cr ;t. Therefore, the joining is only possible between an attractive solution and 
a repulsive one, or between two neutral ones. 

It is easy to see that we can also insert a slice of arbitrary thickness of the 
vacuum solution (|22p between them, obtaining a full relativistic plane "gravitational 
capacitor" . For example, we can trap a slice of Minkowski's space-time between two 
solutions with n = k ct h. 
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7. Attractive Slab surrounded by two different vacuums 

We have already seen that, in the case K CX xt > n > 0, the pressure also vanishes 
inside the slab at the point where u = u . Here we discuss the matching of the slice 
of the interior solution uq < u < n/2 with two vacuums. 
Clearly, the thickness of the slab d is given by 

(tt/2 - u ) 



d = 



(60) 



and • </ . 

Since G(z)\ u=Uo = 1, we can write down from 
in terms of d and p 



the expression which gives k 



^crit 



i + ^d 2 



sin(^/67r pd) 
(2np)\ 4 , 



(61) 

for ^f&Tpd < 1 , (62) 



which clearly shows that k — > K cr it as d — > 0. On the other hand, by using (|96p . we 
can write 



cos l / 3 ( 1 /67rpd) / 3 2 , /- — . 
K= — : V Z , - I 1 - = cos 2 (V67rprf) 



sin(v67rp d) 



7 



(63) 



as d — > \/ 2ip • Thus, as d increases from to J-^-i 



which clearly shows that k 

k monotonically decreases from K cr j t to (see Fig (5]). 

Therefore, the maximum thickness ddec that a solution satisfying the dominant 
energy condition can have, satisfies 



^crit 



3 1/3 ( v / 67rpddcc) - 2Fi( y -\,-\] i;sin 2 (y67rp( 
sin( v / 67rp ddec) 



'dec, 



(64) 



A straightforward numerical computation gives \f&irpd& cc — \.h27AA .... Therefore, 
if < d < ddec, the dominant energy condition is satisfied anywhere. Whereas if 
ddec < d < y^j) there is a region inside the slab where p(z) > p. 

Now, by eliminating k by means of (|6ip . the solution can be parameterized in 
terms of d and p, and (|44|) becomes 



S(z) 



2 Fi(-i,-i; i;sm 2 ( v / 67rp;j 



cos 3 (^/6irp d) sin( v / 67T/9z) 



3 1 / 3 (^67rpz) 



sin(V67rpd) cos3 ( v / 67rp,z) 



2 F 1 (-i,-i;i;sin 2 ( v / 6¥pd)j s i n (76^z) 



sin( v / 67rp d) cos^ (y/Wp z) ^ 

Notice that it clearly shows that (/(— d) = G{0) = 1- By means of (f2"T|) and (|55|) p(z) 



can also be explicitly written down in terms of d and p. The inner line element (|45 
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reads 



-oo < t < oo, 



ds 2 

-oo < X < oo 



Q{zf dt 2 + cos* (V67TP z) (dx 2 +dy 2 ) +dz 2 , 

7T 



-oo < y < oo, 




< -d < z < 0. (66) 



We must impose the continuity of the components of the metric and their first 
derivatives at both matching boundaries, i.e. z — and z = —d. 

The matching at z — was already discussed in section [7] Thus, the upper 
exterior solution, i.e. for z > 0, is the Rindler's space-time 

ds 2 = -(1 + g u z) 2 dt 2 + dx 2 + dy 2 + dz 2 , 



— oo < t < oo, — oo < x < oo, — oo < y < oo, < z < oo , (67) 

which describes a homogeneous gravitational field — g u in the vertical (i.e., z) di- 
rection. And, according to ([51]) . we see that the continuity of d z gu at the upper 
boundary yields 

9u = VSttp (ttcrit - k) , (68) 

which relates the upper external gravitational field g u with matter density p and k. 
By using (|6ip . we can also write it in terms of d and p 

9u= sm(V&Fp-d) ^ F ir\^\'\' Sin2{ ^~ pd) ) - C0S HV^pd)^j . (69) 

At the lower boundary, we have gu{— d) = —Q{—d) 2 = —1, g X x(—d) = g yy (—d) = 
cos 3 (y/Gnpd), and p(— d) = 0. 

On the other hand, regarding the derivatives, since Q{z)\ u=Ua = 1 and 
p{z)\u=u Q = 0, from © we get 

G'(z)\ u=U0 = ~V'{z)\ u=Um = -^JL cotu = tan(V6^c?) , (70) 
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where we have made use of ([T3"]) and Thus, 

S*fftt(-d)|taterior = -2G(-d)G'(-d) = tan(V6^pd) ■ (71) 

While from ([oo]) , we get 

9 z 5 M (-d)|intorior = d z g yy (-d)\ interior = 4 cos 3 ( \/ 6 7175 z) sin(^67T j O (2) .(72) 

Taking into account the discussion in section [4j we can write the corresponding 
lower exterior solution, i.e. for z < —d, as 

ds 2 = - (1 + 3g t (d + z))~% dt 2 + C d (1 + 3gi(d + z))* (dx 2 + dy 2 ) + dz 2 , 

— oo < £ < oo, — oo < x < oo, — oo < y < oo, — d — - — <z<— d, (73) 

which describes a homogeneous gravitational field in the vertical direction and 
finishes up at an empty singular boundary at z — — d — j^. 

Since | exterior = -1 and ^(-d) I exterior = 9yy(-d)\ exterior = C d , we See 

that, taking into account (I45[) . the continuity of the metric components is assured if 
we set Cd = cos 3 (^/6irp d). And, concerning the derivatives of metric's components, 
we have 

d z gtt(z)Ut C ri m = 2gi (l + 3gi{d + z))~ § . (74) 

Therefore, by comparing with (|7Tj) . we see that the continuity of d z gu at the lower 
boundary yields 

gi = tan( v /67rpd) , (75) 

which relates the lower external gravitational field gi with d and p. 
On the other hand, we get from (|T3"|) 

i i - 

^feC^) I exterior = d z g y y(z)\ ex t CT ior = 4 COS 3 ( VOTTp d) #J (1 + 3g/(d + z)) 3 . (76) 

Taking into account ([75]) , by comparing the last equation with ([72]) , we see that the 
matching is complete. 

This solution is remarkably asymmetric, not only because both external gravita- 
tional fields arc different, as we can readily see by comparing ([69]) and ([75]) (see also 
Fig©, but also because the nature of vacuums is completely different: the upper 
one is flat and semi-infinite, whereas the lower one is curved and finishes up down 
bellow at an empty repelling boundary where space-time curvature diverges. 

This exact solutions clearly show how the attraction of distant matter can shrink 
the space-time in such a way that it finishes at an empty singular boundary, as 
pointed out inQ} 

Free particles or photons move in the lower vacuum (— d — ^ < z < — d) along 

the time-like or null geodesies discussed in detail inH. So, all geodesies start and 
finish at the boundary of the slab and have a turning point. Non-vertical geodesies 
reach a turning point point at a finite distance from the singularity (at z — — d— 3^-), 
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and the smaller their horizontal momentum is, the closer they get the singularity. 
The same occurs for vertically moving particles, i. e., the higher the energy, the 
closer they approach the singularity. Only vertical null geodesies just touch the 
singularity and bounce (see Fig. 1 of Ref. ^ upside down) . 

8. The Newtonian limit and the restoration of the mirror 
symmetry 

It should be noted that the solutions so far discussed are mirror-asymmetric. In 
fact, it has been shown in § that the solution cannot have a "plane" of symmetry in 
a region where p(z) > 0. In order to see this, suppose that z = z s is that "plane", 
then it must hold that Q' — V = p 1 = at z s , and so from |9]) we get that also 
p(z s ) = 0, and then Q(z s ) = 1. Now, by differentiating (flQ|) and using (|32|) . we 
obtain p"(z s ) = -Anp 2 < 0. 

Notice that (fT3|) and the condition V'(z s ) = imply that u — ir/2. And then we 
get from (f34|) that the condition G'(z s ) — implies k ~ K cr j t . 

Therefore, the only mirror symmetric solutions is the joining of two identical 
neutral slabs discussed in section [6l We clearly get from (|44]l . that for this solution 
we have 

2Fi(-\,-\; i;sin 2 (V6^z)) 
(l-sin 2 (V6^z)) 1 

which shows that it is a C°° even function of z in —./th- < z < ./t£-. But, of 

course, we have seen in section [2] that — p < p(z) < in this case. 

However, for the solution of the preceding section, this asymmetry turns out to 
disappear when ^/&¥pd <C 1. In fact, from (f69|) we get 

2 



g u = 2npd{l + ^npd 2 + ...) for y^pd < 1 , (78) 

while from (|75|) we get 

gi = 2irpd(l + 2wpd 2 + ...) for ^Tpd < 1 . (79) 

Hence, both gravitational fields tend to the Newtonian result 2wpd, and the differ- 
ence between them is of the order (^/WrFpd) 3 . 
Furthermore, in this limit, (|65|) becomes 

4 JJ 



g(z) - 1 + 2tt P z(z + d) + -'K 1 p 1 z{z i + d A ) + 0((y/6TTpd) b ) , (80) 
so 

gtt(z) = -Q{zf w - (1 + 4npz(z + d)) , (81) 
which shows that the Newtonian potential inside the slab tends to 



$(z) = 2npz(z + d) 



(82) 
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Since $(—5 — z) = $(—5 + z), it is mirror-symmetric at z = —d/2. 

Moreover, we obtain from (fTTj) that the pressure inside the slab tends to the 
hydrostatic Newtonian result 

p(z) = -2np 2 z(z + d) . (83) 



It should also be noted, by comparing ([38)) and (|62|) , that in this limit, they lead to 

^ = ^d 2 . (84) 
P 2 

Therefore, in the Newtonian limit, the the mirror symmetry at the middle point 
of the slab is restored. 



9. Thinner Repelling Slabs 

By exchanging the place of matter and vacuum, we can also match the piece of the 
interior solution discarded in section [7] to the discarded asymptotically flat tail of 

Taub's vacuum, thus getting a repulsive slab. 

Clearly, the inner solution is given by ([66]) . but now — ^/^fp < z — — d. While 
the outer one is given by ([75)) with — d < z. Therefore, we get from ([75)) 



9= ^f^ tan(v^d), (85) 

which relates the external gravitational field g with d and p. But now, the thickness 
of this slab is 



and so 



6lTp cot( v /6^d'). (87) 



3 

Of course, by means of ([27)) , ([65)) and ([86)) . Q{z) and p(z) can also be explicitly 
written down in terms of d' and p in this case. 

In this repulsive case, free particles or photons move in the vacuum (z > 0) along 
the mirror image of time-like or null geodesies discussed in detail in W. All occurs 
in the vacuum as if there were a Taub singularity inside the matter at a distance 
|l/3g| from the surface — this image singularity should not be confused with the 
"real" inner one situated d' from the surface. Therefore, only the Taub's geodesies 
for which the distance between the turning point and the image singularity is smaller 
than should be cut at slab's surface. For instance, this always occurs for vertical 

photons. These facts are easily seen by looking at Fig. 1 of reference ^ upside down 
and by exchanging the position of vacuum and matter. 

Notice that these slabs turn out to be less repulsive than the ones discussed in 
section [5) since all incoming vertical null geodesies reach the slab surface in this 
case. 
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10. Concluding remarks 

We have done a detailed study of the exact solution of Einstein's equations cor- 
responding to a static and plane symmetric distribution of matter with constant 
positive density. By matching this internal solution to vacuum ones, we showed 
that different situations arise depending on the value of a parameter k. 

We found that the dominant energy condition is satisfied only for k > Kdec = 



As a result of the matching, we get very simple complete (matter and vacuum) 
exact solutions presenting some somehow astonishing properties without counter- 
part in Newtonian gravitation: 



The maximum depth that these slabs can reach is y 2lp an< ^ the solutions turn 
out to be remarkably asymmetric. 

We found repulsive slabs in which negative but bounded (\p\ < p) pressure dom- 
inate the attraction of the matter. These solutions finish deep below at a singularity 
where p = —p. If their depth is smaller than y^^f-^, the exterior is the asymptoti- 
cally flat tail of Taub's vacuum plane solution, while when they reach the maximum 
depth the vacuum turns out to be a flat Rindler space-time with event horizons, 
showing that there are incoming vertical photons which never reach the surface of 
the slabs in this case. 

We also found attractive solutions finishing deep below at a singularity. In this 
case the outer solution in this case is a Rindler space-time. 

We also described a non-singular solution of thickness d surrounded by two vac- 
uums. This solution turns out to be attractive and remarkably asymmetric because 
the nature of both vacuums is completely different: the "upper" one is flat and 
semi-infinite, whereas the "lower" one is curved and finishes up down below at an 
empty repelling boundary where space-time curvature diverges. The pressure is pos- 
itive and bounded, presenting a maximum at an asymmetrical position between the 
boundaries. We explicitly wrote down the pressure and the external gravitational 
fields in terms of p and d. We show that if < ^J&npd < 1.52744 . . . , the dominant 
energy condition is satisfied all over the space-time. We also show how the mirror 
symmetry is restored at the Newtonian limit. These exact solutions clearly show 
how the attraction of distant matter can shrink the space-time in such a way that 
it finishes at an empty singular boundary, as pointed out in ^4 

We have also discussed matching an attractive slab to a repulsive one, and 
two neutral ones. We also comment on how to assemble relativistic gravitational 
capacitors consisting of a slice of vacuum trapped between two such slabs. 



0.3513.... 
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Appendix: Some properties of 2 Fi(a, 6; c; x) 

Here, we show how the integral appearing in the first line of (| 16[) is performed. By 
doing the change of variable t = sin 2 u', we can write 

r a i b , , / i r in2 " — * , ,ti i /o+i hi 

/ sin u cos u = - / t 2 (1 - t) 2 dt = -B sin 2 u 
Jo 2 J 2 

where B x (p, q) is the incomplete beta function, which is related to a hypergeometric 
function through 

B x {p,q) = — 2 F 1 (p,l-q;P + l;x) (89) 

(see for example H^j). Therefore, 

' ■ a , bui (sinu) a+1 /a + ll-fca + 3 . 2 
sm a u cos u du = - 2 Fi[ — — , - ; - ; sin u 



o+l V 2 ' 2 ' 2 

= (sin^ 1 (cos^ 1 ^i+^ji + ^jsin^) , (90) 

where we used the transformation 2 .Fi (a, b; c; x) = (1 — x) c ~ a ~ b 2 F 1 (c — a,c — b; c; x) 
in the last step. 

For the sake of completeness, we display here the very few formulas involving 
hypergeometric functions 2 Fi(a, b;c;z) required to follow through all the steps of 
this paper. 

As it is well known 

t-, / , n ., a b a(a + 1)6(6 + 1) z 2 . 
2 F 1 (a,b;c;z) = l + —z+ c(c+1) -^[ + --- . for 1*1 < 1 ■ ( 91 ) 

By using the transformation H 2 | 13 | 

2 F 1 (a,6;c;z) = ) C> {C ~ ° - 2^(0, 6; o + 6 - c+ 1; 1 - z) 
1 (c — a)I (c — 0) 

+ {1 _ z) c-u-b r(c)r(o + 6-c) 6;c _ a _ 6 + 1 (92) 

r(a)r(6) 

with a = —1/2, 6 = —1/6, and c = 1/2, we find the useful relations 

._^|)^ +1+ | + |! + il|! + ..., totW<1 , (M) 

or, by making z — > 1 — z, 
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